Reducing or enhancing chaos using periodic orbits 
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A method to reduce or enhance chaos in Hamiltonian flows with two degrees of freedom is dis- 
cussed. This method is based on finding a suitable perturbation of the system such that the stability 
of a set of periodic orbits changes (local bifurcations). Depending on the values of the residues, 
reflecting their linear stability properties, a set of invariant tori is destroyed or created in the neigh- 
borhood of the chosen periodic orbits. An application on a paradigmatic system, a forced pendulum, 
illustrates the method. 

PACS numbers: 05.45.-a 



Changing the dynamical properties of a sys- 
tem is central to the design and performance of 
advanced devices based on many interacting par- 
ticles. For instance, in particle accelerators, the 
aim is to find the appropriate magnetic elements 
to obtain an optimal aperture in order to increase 
the luminosity of the beam, thus requiring the 
decrease of the size of chaotic regions. In plasma 
physics, the situation is slightly more complex : 
Inside a fusion device (like a tokamak or a stel- 
larator), one needs magnetic surfaces in order to 
increase confinement. These surfaces are invari- 
ant tori of some fictitious time dynamics. A con- 
trol strategy would be to recreate such magnetic 
surfaces by an appropriate modification of the ap- 
paratus (magnetic perturbation caused by a set of 
external coils). On the opposite, in order to col- 
lect energy and to protect the wall components, 
an external modification of the magnetic equilib- 
rium has to be performed such that there is a 
highly chaotic layer at the border (like an ergodic 
divertor). Therefore these devices require a spe- 
cific monitoring of the volume of bounded mag- 
netic field lines. Another example is afforded by 
chaotic advection in hydrodynamics : In the long 
run to achieve high mixing in microfluidics and 
microchannel devices in particular, the presence 
of regular region prevents such mixing, and hence 
a possible way to enhance mixing is to perturb ex- 
ternally the system according to some theoretical 
prescriptions, in order to destroy invariant sur- 
faces. 
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I. INTRODUCTION 

A very fruitful information on the dynamics can be 
gained from the study of periodic orbits pj. First, be- 
cause these particular orbits are generically almost ev- 
erywhere in phase space, and second because they can 
be computed easily, i.e. with some short integration time. 
These periodic orbits together with their stability orga- 
nize locally the dynamics. It is then natural to consider 
them as a cornerstone of control strategies. For instance, 
in order to create invariant tori of Hamiltonian systems, 
Cary and Hanson 0, Q proposed a method based on the 
computation of an indicator of the linear stability of a 
set of periodic orbits, namely Greene's residue p|. It 
provides an algorithm to find the appropriate values of 
some pre-defined parameters in order to reconstruct in- 
variant tori by vanishing some selected residues. First de- 
veloped for two-dimensional symplectic maps, it has been 
extended to four dimensional symplectic maps, and has 
been applied to stellarators [j| (where periodic orbits are 
closed magnetic field lines) and particle accelerators @- 

In this article, we review and extend this residue 
method. The aim is to tune appropriately the param- 
eters of the system such that appropriate bifurcations 
occur. It is well-known in the literature that local bifur- 
cations occur when the tangent map associated with the 
Poincare map obtained by a transversal intersection of 
the flow, has an eigenvalue which is a root of the unity. In 
particular, periodic orbits can lose their stability in case 
of multiple eigenvalues on the unit circle, i.e. when these 
eigenvalues are equal to 1 or — 1 for two-dimensional sym- 
plectic maps. Therefore it is natural to consider Greene's 
residues as a way to locate those bifurcations. In this con- 
text, vanishing residues indicate the specific values of the 
parameters where significant change occurs in the system 
and hence will be the basis for the reduction of chaos (by 
creation of invariant tori) as in Refs. 0, but also for 
the destruction of regular structures. 

In Sec. [n] we review some basic notions on periodic 
orbits of Hamiltonian systems and their stability, and we 
explain the details of the residue method. We give the 
condition on the residues of a pair of Birkhoff periodic 
orbits to create an invariant torus in their vicinity, and a 
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similar condition which leads to a destruction of nearby 
invariant tori. In Sec. IIIII we apply this method to the 
destruction and creation of librational and rotational in- 
variant tori of a particular Hamiltonian system, a forced 
pendulum with two interacting primary resonances, used 
as a paradigm for the transition to Hamiltonian chaos. 

II. THE RESIDUE METHOD 

We consider an autonomous Hamiltonian flow with two 
degrees of freedom which depends on a set of parameters 
denoted a — (ax, 02, . . . , a m ) € M m : 

i = IVH(z;a), 

where z = (p, q) S M. 4 and J = ( ® J 2 ~\ , and I2 

being the two-dimensional identity matrix. In order to 
determine the periodic orbits of this flow and their lin- 
ear stability properties, we also consider the tangent flow 
written as 

^-J t (z) = $V 2 H(z;a)J t , 

where J° = I 4 and V 2 H is the Hessian matrix (composed 
by second derivatives of H with respect to its canonical 
variables). For a given periodic orbit with period T, the 
spectrum of the monodromy matrix J T gives its linear 
stability property. As the flow is volume preserving, the 
determinant of such a matrix is equal to 1. Moreover, if 
A is an eigenvalue, so are 1/A, A* and 1/A*. As the orbit 
is periodic, A = 1 is an eigenvalue with an eigenvector 
in the direction of the flow. Its associated eigenspace is 
at least of dimension 2 since there is another eigenvector 
with eigenvalue 1 coming from the conserved quantity 
H = E. Therefore, according to the remark above, the 
orbit is elliptic if the spectrum of J T is (1, 1, e M , e~ luJ ) 
(and stable, except at some particular values), or hyper- 
bolic if the spectrum is (1,1, A, 1/A) with A G M* (un- 
stable). The intermediate case is when the spectrum is 
restricted to 1 or —1 and the orbit is called parabolic. 
Whether or not the parabolic periodic orbit is stable de- 
pends on higher order terms. In a more concise form, the 
above cases can be summarized using Greene's definition 
of a residue which led to a criterion on the existence of 
invariant tori 0,0 : 



We notice that the 4 (instead of 2 for 2D maps) in the 
numerator comes from the two additional eigenvalues 1 
coming from autonomous Hamiltonian flows. If R <e]0, 1[, 
the periodic orbit is elliptic; if R < or R > 1 it is 
hyperbolic; and if R = and R = 1, it is parabolic 
and higher order expansions give the stability of such 
periodic orbits. Since the periodic orbit and its stability 
depend on the set of parameters a, the features of the 



dynamics will change with variations of the parameters. 
Generically, periodic orbits and their linear stability are 
robust to small changes of parameters, except at specific 
values where bifurcations occur. The proposed residue 
method to control chaos detects these rare events to yield 
the appropriate values of the parameters leading to the 
prescribed behavior on the dynamics. 

The residue method which leads to a reduction or an 
enhancement of the chaotic properties of the system is 
based on the change of stability of periodic orbits upon 
a change of the parameters of the system. For a = 0, 
let us consider two associated Birkhoff periodic orbits 
(i.e. periodic orbits having the same action but different 
angles in the integrable case and having the same rotation 
number on a selected Poincare section), one elliptic O e 
and one hyperbolic Oh- Let us call R e and Rh their 
residues. We have R e (0) > (and smaller than one) 
and Rh(0) < 0. We slightly modify the parameters a 
until the elliptic periodic orbits becomes parabolic. Some 
particular situations arise at some critical value of the 
parameters a. = a c : 

(l) : R e (a c ) = R h (a c ) = 0. 

(ii) : R e (a c ) = while Rh(a c ) < 0. 
(Hi) : R e (a c ) = 1 while Rh(a c ) < 0. 

The first case is associated with the creation of an in- 
variant torus. The two latter cases might be associated 
with the destruction of invariant tori (the ones around 
the elliptic periodic orbit). The third one is associated 
with a period doubling bifurcation. In this latter case, 
the change of stability of the new elliptic periodic or- 
bit has to be considered. Other interesting cases occur 
depending on the set of selected periodic orbits. The 
situation (i) resembles the integrable situation where all 
the residues of periodic orbits of constant action are zero. 
It is expected that an invariant torus is reconstructed in 
this case. It can be associated with a transcritical bifur- 
cation (an exchange of stability) , a fold, or another type 
of bifurcation. In the situation (ii), a change of stabil- 
ity occurs : The elliptic periodic orbit turns hyperbolic 
while the hyperbolic one stays hyperbolic. It is generi- 
cally characterized by a stationary bifurcation. In this 
case, the destruction of invariant curves is expected in 
general whether there are librational ones (representing 
the linear stability of an elliptic periodic orbit) or the 
neighboring rotational ones. 

An extra caution has to be formulated since this 
method only provides an indicator of the linear stability 
of periodic orbits. The nonlinear stability (or instabil- 
ity) has to be checked a posteriori by a Poincare section 
for instance. This method only states that a bifurcation 
has occurred in the system, whether it is a stationary, 
transcritical, period doubling or other types of bifurca- 
tions. A more rigorous and safer control method would 
require to consider the global bifurcations, like the ones 
obtained by the intersections of the stable and unstable 
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manifolds of two hyperbolic periodic orbits in the spirit 
of Ref. @. However such a control method would be 
computer-time consuming (determination of the stable 
and unstable manifolds) and hence not practical if some 
short time delay feedback is involved in the control pro- 
cess. 



III. APPLICATION TO A PARADIGMATIC 
MODEL 

We consider the following forced pendulum system 
with 1.5 degrees of freedom 

P 2 

H(p, x, t) = — + e (cosx + cos{x — t)) . (1) 

A Poincare section of Hamiltonian (JJJ is depicted on 
Fig.[T]for s = 0.065 and on Fig.Hfor e = 0.034. In order 
to modify the dynamics of Hamiltonian l|T]l. we add an 
additional (control) parameter k : We consider a family 
of Hamiltonians of the form 

H c (p, x,t) = he (cos a; + cos(x — t))-\ — e 2 cos(2x — i), 

where k is not too large in order to consider a small mod- 
ification of the original system, and minimizing the en- 
ergy cost needed to modify the dynamics. Other choices 
of families of control terms are possible (not restricted to 
k cos(2x~t)). In particular, more suitable choices of con- 
trol terms would include more Fourier modes. We have 
selected a one-parameter family which originates from 
another control strategy which has been proved to be ef- 
fective 0- The goal here is to determine the particular 
values of the parameter k such that suitable modifications 
of the dynamics (which will be specified later) occur. 

The algorithm is as follows : First, we determine 
two periodic orbits of Hamiltonian Q, an elliptic and 
a hyperbolic one with the same rotation number on 
the Poincare section, using a multi-shooting Newton- 
Raphson method for flows [ij. Then we modify contin- 
uously the control parameter k and follow these two pe- 
riodic orbits. We compute their residues as function of 
k. 

A first analysis is done on librational invariant tori 
(around the primary resonance located around p rs 0). 
We point out in Fig. ^ three particular elliptic periodic 
orbits (and their associated hyperbolic ones), labeled 1,2 
and 3, with, respectively, Q — 4, Q = 9 and Q = 13 
intersections with the Poincare section (t = mod 2tt). 
These orbits will be used for two purposes : First we 
follow the idea of Cary and Hanson on the construction 
of invariant tori. Then we extend the residue method to 
the destruction of these tori. A similar analysis is done 
on rotational invariant tori (the example of the golden- 
mean invariant torus is treated). This case allows us to 
compare the residue method with another approach on 
the control of Hamiltonian systems. 



1.5 




-0.5\ | | | ■ | 

1 2 3 4 5 6 

X 



a. 







■ 2 

X x 




—» pet • 


3 

X 










* X " 


• • 






" X x X • 





1 2 3 4 5 6 



X 

FIG. 1: Poincare sections of Hamiltonian Q with e = 0.065. 
The arrows indicate the elliptic periodic orbits for the three 
cases considered here. 
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FIG. 2: Poincare section of Hamiltonian Q for e = 0.034. 



Briefly we determine the control parameter k such that 
there is a creation of an invariant torus if the original 
system does not have one, and the destruction of an in- 
variant torus if the system does have one. 
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FIG. 3: Residues R of the elliptic and hyperbolic (bold line) 
periodic orbits of Case 1 (Q — 4) as functions of the parameter 
k for Hamiltonian with e = 0.065. 



A. An exchange of stability associated with the 
creation of librational invariant tori 

For the case Q = 4, Fig. [3] represents the values of 
the residues of the elliptic and hyperbolic periodic or- 
bits as functions of the control parameter k [see Eq. PJl]. 
At k — k c w 2.747, both residues vanish which means 
that they become parabolic periodic orbits as in the in- 
tegrable case. By increasing k, we notice that both orbits 
exchange their stability which is the manifestation of a 
transcritical bifurcation while each of the periodic orbits 
undergo individually a tangent bifurcation. This type 

of bifurcation has been observed in Rcfs. nam. At 

k = k c , an invariant torus is reconstructed. In order to 
check the robustness of the method, one could argue that 
since this invariant torus is composed of periodic orbits, 
it is not expected to be robust. However, by continuity 
in phase space, an infinite set of invariant tori is present 
in the neighborhood of the created invariant torus. Most 
of them have a frequency which satisfy a Diophantine 
condition and hence which will persist under suitable hy- 
pothesis on the type of perturbations. 

The locations of the different periodic points on the 
Poincare section as k varies are indicated by arrows on 
Fig. 0] The change of stability of these periodic points 
is associated with the creation of an invariant torus (also 
represented in Fig.0|by the plot of the separatrices). We 
notice that apart from the exchange of stability, the phase 
space in the neighborhood of these periodic orbits is still 
regular (the chaotic region around the hyperbolic peri- 
odic orbits is not well developed), and hence the regular 
nature of phase space has not been changed locally (or 
one needs to consider higher values of the parameters) . 

The same analysis can be carried out on a set of pe- 
riodic orbits which are located in a more chaotic region, 
like for instance the two cases Q = 9 and Q = 13 of 
Fig. outside the regular resonant island. The values 
of the residues as functions of the parameter k are re- 
spectively represented on Figs. [S] and for Q — 9 and 
Q = 13. For Q — 9, we notice that the residues do not 
vanish in the range of k we considered although there are 
small and extremum at the same value of the parameter 
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FIG. 4: Poincare section of periodic orbits and some trajecto- 
ries for Case 1 (Q — 4) and Hamiltonian J2J with e = 0.065. 
The trajectories in gray are for k — and the ones in black 
are for k = 5. At k = kc, an invariant torus of the system 
is represented (bold line). The arrows indicate the change of 
locations of the periodic points as k increases. 
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FIG. 5: Residues R of the elliptic and hyperbolic periodic 
orbits of Case 2 (Q = 9) as functions of the parameter k for 
Hamiltonian with e = 0.065. 

k as 2.79. However, even if these residues do not vanish 
(and therefore no exchange of stability by the creation of 
an invariant torus), there is a significant regularization 
of the dynamics at this specific value of the parameter 
(not shown here). For Q = 13, the residues vanish for 
k w 2.76 (see Fig.[f)J) and there is a transcritical bifurca- 
tion associated with the creation of a set of invariant tori 
like in Figs. and 0] The associated phase space shows 
a significant increase of the size of the resonant island 

B. Enhancing chaos near a resonant island 

In this section, we address the destruction of a reso- 
nant island by breaking up librational invariant tori. We 
notice that on Fig. EI a bifurcation occurs at k w —1.254 
for the Case 2 (Q = 9) when the residue of the elliptic pe- 
riodic orbit becomes equal to 1. The neighborhood of this 
periodic orbit becomes a chaotic layer and the qualita- 
tive change in the dynamics is seen since the chaotic layer 
becomes thicker at this value of the parameter. However 
since this periodic orbit was initially (at k — 0) already in 
the outer chaotic region (see Fig.^), the regularization is 
not drastic. In order to obtain a more significant change 
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FIG. 6: Residues R of the elliptic and hyperbolic periodic 
orbits of Case 3 (Q = 13) as functions of the parameter k for 
Hamiltonian @ with e = 0.065. 

in the dynamics and a large chaotic zone, one needs to 
select a periodic orbit inside a regular region, like for in- 
stance the one with Q — 13. A bifurcation occurs at 
fe w —2.484 where the residue of the elliptic periodic or- 
bit crosses 1 (see Fig-EJ ■ A Poincare section for the latter 
case is depicted on Fig. and shows that a significantly 
large neighborhood has been destabilized by the control 
term. We notice that the ratio between the size of the 
control term and the one of the perturbation is equal to 
fee w 0.16. We also notice that this last value of fc is 
larger than the one required for Q = 9. As expected, one 
needs a larger amplitude to destabilize a region closer to 
a regular one. A more effective destabilization procedure 
can be obtained with the periodic orbit Q — 4 which is 
inside the regular region. However, as mentioned, the 
value necessary for this destabilization (fe ps —12.5 or for 
fe i=s 16.3) is too large; hence we discard it because of our 
restriction on energy cost. 



C. Creation of the goldenmean rotational invariant 

torus 

In this section, we apply the same approach on rota- 
tional invariant tori. It allows us to compare the results 
with the ones obtained by a control method proposed in 
Refs. 0,0,0- First, the idea is to look at the creation 
of a specific invariant torus. For instance, we select a 
torus which has been widely discussed in the literature 
(see Ref. ^4| and references therein), the goldenmean 
one, which has a frequency w = (3 — y/E) /2 for Hamilto- 
nian We choose e = 0.034, and we first notice that 
when k — 0, this Hamiltonian does not have such an in- 
variant torus (since its critical value is e ~ 0.02759 [l4T|). 
The purpose here is to find the value of the control pa- 
rameter k needed by the residue method to reconstruct 
this invariant torus (such that Hamiltonian (J2J) has this 
invariant torus). 

The idea of doing this follows Greene's residue crite- 
rion. By performing an appropriate change of stability 
on higher and higher order periodic orbits, the ampli- 
tude of the control term should be smaller and smaller. 
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FIG. 7: Poincare section of Hamiltonian for e = 0.065 and 
k = -2.484. 

For Q = 2, the residues of the elliptic and hyperbolic 
periodic orbit vanish at fc « 4.3, and for Q = 3, at 
k w 3.4. For Q = 5, both residues vanish at k = 2.98 
and also at fc « 5.07. At these values of the parameter, 
the phase space is locally filled by invariant tori where it 
is also expected that the goldenmean invariant torus is 
present. We notice that the elliptic periodic orbit with 
Q = 8 (the next one in Greene's residue approach for 
the analysis of the golden mean torus) is destabilized at 
e pa 0.0325. Therefore, there is no elliptic periodic orbit 
with Q — 8 at e — 0.034 and the analysis using the cou- 
pled elliptic/hyperbolic periodic orbits cannot be carried 
out. However, by following the two (initially hyperbolic) 
periodic orbits with Q = 8, we see that both residues 
vanish at k 2.84. 

We compared these values of stabilization with the one 
given by a method of local control based on an appropri- 
ate modification of the potential to reconstruct a specific 
invariant torus 0, Il2| . Such method provides explic- 
itly the shape (and amplitude) of possible control terms 
whereas the one used in this article has been guessed from 
these references. By appropriate truncation (keeping the 
main Fourier mode), this method provides 

£ 2 

f( x > *) = 771 \ cos ( 2x - *)> 

2w(l — ll>) 

where uj — (3 — \/5)/2, as an approximate control term. 
Therefore the amplitude is fe = l/w(l — w 4.24 which 
is of the same order as the values obtained by zeroing the 
residues. However, we point out that smaller values are 
obtained by looking at higher periodic orbits. 

Therefore, an efficient control strategy is to combine 
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the advantages of both methods : First, the specific shape 
of the terms that have to be added to regularize the sys- 
tem is obtained using the method of Ref . ■ Then the 
amplitudes of these terms are lowered using high order 
periodic orbits. By considering the control term used in 
this article, we expect that zeroing the residues of high 
period will not be feasible with just this term (as it is 
the case for instance in Fig.[5J). A more suitable form of 
control terms would be constructed from an exact control 
term which is 

= mi-") 008(28 "*> 

e 2 e 2 

-—r cos 2a; - — ^ cos2(ir - t). 

Au 2 4(1 -lu) 2 y ' 

However, it should be noticed that a control term given 
by Ref. [|| is not always experimentally accessible. The 
idea is to use a projection of this control term onto a 
basis of accessible functions. This projected control term 
would give an idea of the type of control terms to be used 
for the residue method. 

We would like to stress that in the absence of elliptic is- 
lands an initial guess for the Newton-Raphson method is 
not straightforward from the inspection of the Poincare 
section. In particular, it is not easy to select the ap- 
propriate hyperbolic periodic orbits which will lead to a 
significant change in the dynamics. However, once it has 
been located, the method can follow them by continuity 
in the same way as the elliptic ones since the Newton- 
Raphson method does not depend on the linear stability 
of these orbits. This makes the method more difficult 
(although possible) to handle for just hyperbolic periodic 
orbits. 



D. Destruction of the goldenmean rotational 
invariant torus 

In this section, we consider Hamiltonian J2J with e = 
0.0275. We notice that for k = 0, Hamiltonian (0) docs 
have the rotational goldenmean torus. The purpose is 
to find some small values of the parameter k where this 
invariant torus is destroyed. We notice that this case is 
easier to find than in the previous section since it is well- 
known that any additional perturbation will end up by 
destroying an invariant torus generically. Here it means 
that there will be large intervals of parameters for which 
the torus is broken (contrary to the case of the creation 
of invariant tori). However we will add an additional 
assumption that the parameters for which this invariant 
torus is destroyed has to be small compared with the 
perturbation. We also notice that the destruction of the 
golden mean invariant torus is first obtained for negative 
values of the control parameter (see Fig- - 

First we illustrate the method by considering specific 
elliptic and hyperbolic periodic orbits (with winding ratio 
5/13) near the goldenmean torus which will show the 
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FIG. 8: Residues R of the elliptic and hyperbolic periodic 
orbits with Q — 13 and also of the one with Q — 26 (dashed 
line) born out of a period doubling bifurcation for Hamilto- 
nian © with e = 0.0275. 



changes of dynamics occurring as the parameter is varied. 
We notice that the behaviors described below are generic 
for all the neighboring periodic orbits. 

The residues of these periodic orbits as functions of the 
parameter k are shown in Fig. |SJ We notice that the el- 
liptic periodic orbit changes its stability, i.e. becomes hy- 
perbolic, at k w —1.205 (where its residue becomes equal 
to 1). A close inspection of the Poincare section shows 
on Fig. Ej] that it undergoes a period doubling bifurcation 
into an elliptic periodic orbit with 26 intersections on the 
Poincare section (and winding ratio 10/26) which has a 
residue zero at the bifurcation. By following the residue 
of this elliptic periodic orbit (depicted by a dashed line 
in Fig. [HJ we see that it vanishes for k = —1.6256. At 
this value of the parameter and for higher value in am- 
plitude, all the periodic orbits considered here (the two 
with Q = 13 and the one with Q = 26) are hyperbolic. 
Therefore it is expected that there is a chaotic zone in this 
area and it is a value at which the torus is expected to be 
broken (confirmed by a close inspection of the Poincare 
section). 

It is important to notice that a vanishing residue does 
not automatically imply that there is a creation of an in- 
variant torus, contrary to the previous cases which were 
obtained by using jointly the elliptic and hyperbolic pe- 
riodic orbits (and vanishing residues in both cases) . Here 
the hyperbolic periodic orbit associated with these ellip- 
tic periodic orbits (which is the periodic orbit from which 
the new elliptic orbit was born out by a period doubling 
bifurcation) stays hyperbolic as the residue of the elliptic 
one vanishes. This feature is generic : The same analysis 
has been carried out for higher order elliptic periodic or- 
bits close to the goldenmean invariant torus, i.e. the ones 
with winding ratio 8/21, 13/34, 21/55, 34/89 : First, the 
values of the control parameter for which the residues 
(which are around 0.25 for k = and increase as k de- 
creases) cross 1 are computed and reported in Table I 
(denoted k(R — 1)). At these values of the parameters, 
a period doubling bifurcation occurs for each of them. 
Then we follow the residues of the elliptic periodic or- 
bits with double period Q — 42, Q = 68, Q — 110 and 
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FIG. 9: Poincare section around the periodic orbit with wind- 
ing ratio 5/13 (indicated with crosses) for Hamiltonian for 
e = 0.0275 and k = —1.215. The period orbit with period 26 
(indicated by circles) results from a period doubling bifurca- 
tion of the one with period 13 (represented by crosses on the 
Poincare section). 



Q 


13 


21 


34 


55 


89 


k{R=l) 


-1.205 


-0.705 


-0.435 


-0.273 


-0.179 


k c 


-1.626 


-0.935 


-0.566 


-0.350 


-0.225 



TABLE I: Values of the parameter k at which the residue of 
the elliptic periodic orbit with period Q crosses 1 (denoted 
k(R — 1)) and at which the residue of the elliptic periodic 
orbit with period 2Q obtained by period doubling bifurcation 
at k(R = 1) vanishes (denoted k c ). 



that the goldenmean invariant torus is destroyed by this 
additional perturbation but not the ones in the neighbor- 
hood. If one is looking at large scale transport properties, 
these other invariant tori have to be taken into account. 



Concluding remarks 

In this article, we reviewed and extended a method of 
control of Hamiltonian systems based on linear stabil- 
ity analysis of periodic orbits. We have shown that by 
varying the parameters such that the residues of selected 
periodic orbits cross or 1, some important bifurcations 
happen in the system. These bifurcations can lead to the 
creation or the destruction of invariant tori, depending 
on the situation at hand. Therefore we have proposed 
a possible extension of the residue method to the case 
of increasing chaos locally. Moreover, we have compared 
two methods of chaos reduction, and by taking advan- 
tage of both methods, we have devised a more effective 
control strategy. It is worth noticing that the extension 
of Cary-Hanson's method to four dimensional symplectic 
maps has been done in Refs. for the increase of dy- 

namic aperture in accelerator lattices. The extension to 
the destruction of invariant surface would be to consider 
the change of linear stability of selected periodic orbits. 
However, it would require to consider new types of bi- 
furcations which occurs in the system, like for instance, 
Krein collisions [lq . 



Q = 178. The parameter values at which these residues 
vanish are also reported in Table I. For instance, using 
the periodic orbit with winding ratio 8/21, we obtain 
k = —0.935 as the value at which the residue of the bi- 
furcated elliptic periodic orbit with winding ratio 16/42. 
If we consider higher order periodic orbits, it happens 
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